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AXISYMMETRIC FLOW OF IDEAL FLUID MOVING IN A 
NARROW DOMAIN: A STUDY OF THE AXISYMMETRIC 
HYDROSTATIC EULER EQUATIONS 


ROBERT M. STRAIN AND TAK KWONG WONG 


Abstract. In this article we will introduce a new model to describe the 
leading order behavior of an ideal and axisymmetric fluid moving in a very 
narrow domain. After providing a formal derivation of the model, we will 
prove the well-posedness and provide a rigorous mathematical justification for 
the formal derivation under a new sign condition. Finally, a blowup result 
regarding this model will be discussed as well. 
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1. Introduction 

In various applications in meteorology, oceanography, atmospheric dynamics, 
blood flow and pipeline transport, the vertical or radial length scale of the 
underlying flow is usually small compared to the horizontal length scale. To study 
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these problems, the standard approach is to apply the hydrostatic approximations. 
For example, when a two-dimensional ideal fluid moves in a fixed and very narrow 
channel, one can describe the leading order behavior of the fluid motion by the 
two-dimensional hydrostatic Euler equations, which can be formally derived by 
the hydrostatic limit or the least action principle [3]. Under the local Rayleigh 
condition, the formal derivation of the two-dimensional hydrostatic Euler equations 
via the hydrostatic limit was rigorously justified in [SlEllig. Without the local 
Rayleigh condition, the formal derivation may not be valid [5JIS]. The local-in-time 
existence and uniqueness under the analyticity assumption [5], the local Rayleigh 
condition cma, or their combinations in different regions [7] are also known, but 
the global-in-time existence is still open. Furthermore, for a general initial data, 
the two-dimensional hydrostatic Euler equations are somewhat ill-posed: see [M] 
for the linearized instability, and min] for the formation of singularities. 

In this paper, we study the leading order behavior of axisymmetric and ideal 
flows moving in a very narrow domain in three spatial dimensions. The prime 
objectives of this paper are as follows: 


(i) to formally derive the axisymmetric hydrostatic Euler equations, which 
describe the leading order behavior of axisymmetric Euler flows moving in 
a thin tube, via the hydrostatic limit (see Subsection mi); 

(ii) to introduce a new sign condition (see inequality (2.8) below), which is an 
analogue of the local Rayleigh condition in two spatial dimensions, for the 
axisymmetric hydrostatic Euler equations in three spatial dimensions; 

(iii) to prove the well-posedness of the axisymmetric hydrostatic Euler equations 
under the new sign condition (see Theorem 12.21 Sections [3] and 01 as well as 
Appendix ICll: 

(iv) to provide a rigorous mathematical justification of the formal derivation for 
the axisymmetric hydrostatic Euler equations under the new sign condition 
(see Theorem 12.41 and Section [5]); 

(v) to discuss the finite time blowup of smooth solutions of the axisymmetric 
hydrostatic Euler equations (see Theorem 12.61 and Section [6]). 


To the best of our knowledge, these issues have not been studied in the literature. 

The main difficulty for the study of axisymmetric hydrostatic Euler equations 
is the loss of the horizontal regularity. Similar to the two-dimensional hydrostatic 
Euler equations, the axisymmetric hydrostatic Euler equations are also derived by a 
singular limit process, called the hydrostatic limit in the literature. The hydrostatic 
limit/approximation usually simplifies the pressure term, but creates a loss of one 
horizontal derivative. Due to the regularity loss in the horizontal direction, standard 
energy methods do not apply in general. 

The main novelty of this work is to introduce a new sign condition, which is an 
analogue of the local Rayleigh condition in two spatial dimensions. Under this sign 
condition, the horizontal regularity loss can be avoided by the nonlinear cancelation 
(2.9) below. As a result, the theory for the axisymmetric hydrostatic Euler 


equations can be established under this new sign condition by using the standard 
energy method. Furthermore, in order to simplify the computations, we will make 
use of new differential operator and dependent variables, which will be stated in 
(3.5) below. Under these new differential operator and dependent variables, the 


vorticity system for the axisymmetric hydrostatic Euler equations is equivalent 
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to that for the two-dimensional hydrostatic Euler equations in a certain sense. 
Therefore, the analysis in this paper is similar to that in mm- 

The rest of this paper is organized as follows. First of all, we will provide the 
formal derivation of the axisymmetric hydrostatic Euler equations, introduce the 
new sign condition (i.e., inequality (2.8)) and state our main results (i.e.. Theorems 
12.2112.41 and EH) in Section El In Section El we will derive a priori estimates, and 
apply these estimates to prove the uniqueness and stability. The existence will 
be shown in Section E) Using the entropy method, we will provide a rigorous 
mathematical justification of the formal derivation for the axisymmetric hydrostatic 
Euler equations in Section [5j Finally, the blowup result will be discussed in Section 
[6l For the sake of self-containedness, we will also provide elementary proofs in 
Appendices lAllCl 

Let us end this introduction by commenting on our notation. Throughout this 
paper, all constants with or without subscript(s) may be different in different lines. 
Unless mentioned otherwise, a constant with subscript (s) illustrates the dependence 
of the constant, for example, Cs^a is a constant depending on s and a only. 


2. Formal Derivation and Main Results 

In this section, we will first provide a heuristic derivation of the axisymmetric 
hydrostatic Euler equations in Subsection 12.11 and then state the main results of 
this paper in Subsection 12.21 


2.1. Formal Derivation. The aim of this subsection is to formally derive the 
axisymmetric hydrostatic Euler equations via a rescaling limit. This rescaling limit, 
called the hydrostatic limit, can be described as follows. 

An ideal fluid moving in a periodic and narrow domairQ := {(A, Y, Z); X G 
T := R/Z, Y'^ + Z'^ < e^} is governed by the usual three-dimensional incompressible 
Euler equations: 

' dtUe + Ue- VU, := -VP, in (0, T) x U, 

div U, = 0 in (0, T) x U, 

(2-1) 

Ue ■ n\dn^ = 0 
Ue\t=0 = UeO, 

where U, is an unknown velocity field, P, is an unknown scalar pressure, Ueo is a 
given initial velocity field, and n is the outward unit normal vector to the boundary 
90 ,. 

Now, let us assume that the underlying flow is axisymmetric without swirl. In 
other words, using the cylindrical coordinates 

X := X, Y := Pcos0, and Z := Psin0, 


we can express the velocity field P, and the pressure P, as follows: 


P, := U^{t, X, R)ex + U^{t, A, R)eR 
P, :=P,(t,A,P), 


^In order to simplify our presentation, we only consider a periodic tube (i.e., X GT := R/Z), 
but one could consider an infinite tube (i.e., X G R) or other physical domains. 
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where ex and cr are the unit vectors in the X and R directions respectively. Under 
the assumption (2.2)[ the usual incompressible Euler equations (2.1) become the 
axisymmetric Euler equations: for {t,X,R) G (0,T) x T x (0,e), 


(2.3) 


dtU^ + U^dxU^ + U^OrU^ = -dxPe 
dtUP + U^dxU^ + U^BrU^ = -OrP, 

dxU^ + ^dR{RU^^) = 0 

K 


u 


\R=e 


= 0 


U^\t^o = U^o- 

Here, we only have to impose the initial condition for because that for can 


be uniquely determined by using (2.3)(2.3) 


In order to study the leading order behavior of the flow as the thickness e goes 
to 0’’’, the standard approach is to rescale the physical domains into a uniform 
domain. More precisely, we apply the rescaling 


(2.4) 


X := X, 

U^{t,X,R) ■= [t,X, 

Pe{t,X,R) := pe (t,X, 


_ R 
e ’ 


U^it,X,R) ■= e< t,X 


R 


{t, X, r) G (0, T) X T X (0,1), 


to rewrite the system (2.3) as the axisymmetric rescaled Euler equations: for 


(2.5) 


= -dxPe 

HdtU: + U^d^ul + uldrUl) = -drPe 


+ -dr{rul) = 0 


<|r=l = 0 

<|i=o = <0> 


where (u^, u^) and p^ are unknowns, and u^q is the given initial horizontal velocity. 

Eormally, if {u^,ul,pe) converges to {u^,u’",p) as e goes to 0'^, then {u^,u'~,p) 
will satisfy the axisymmetric hydrostatic Euler equations: for {t,x,r) G (0,T) x 
Tx (0,1), 

dtu^ + u^dxU^ + u^drU^ = —dxP 


( 2 . 6 ) 


drP — 0 

dxU^ +-driru^) = 0 
r 

vJ'\r=\ = 0 

u^\t=o=ul. 


It is worth noting that equation (2.6), is equivalent to the fact that the scalar 
pressurep is independent of r, i.e., p := p(t, x). Therefore, system (2.6) is equivalent 


to system (2.7) below provided that p is assumed to be independent of r. 
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The system (2.6) formally describes the leading order behavior of axisymmetric 
and ideal flows moving in the narrow domain The above limiting process is 
called the hydrostatic limit. 

Let us end this subsection by explaining why the vorticity w for the axisymmetric 
hydrostatic Euler equations (2.6) (or equivalently, (2.7) I is drU^ as follows: 


Remark 2.1 (Vorticity for the Axisymmetric Hydrostatic Euler Equations). It is 
well-known that the vorticity for the axisymmetric (without swirl) velocity field 
(2.2) ^ is just —{dnU^ — dxU^)ee, where ee is the unit vector in the 0 direction. 


Applying the rescaling (2.4) we have 


- 


dxU,^ = - [druf - e^d^u:) . 


As e —>■ 0+, the quantity druf — e^dxul converges to drU^ formally. Therefore, 
we denote the vorticity for the axisymmetric hydrostatic Euler equations (2.6), 
(or equivalently, (2.7) I as the scalar quantity drU^, which is corresponding to the 
leading order term of the vorticity as e —>■ O'*". 

2.2. Main Results. In this subsection we will introduce the axisymmetric hydro¬ 
static Euler equations, the new sign condition and our main results. 

According to the formal derivation in Subsection 12.11 the leading order behavior 
of an axisymmetric (without swirl) ideal flow moving in a periodic and narrow 
channel can be described by the axisymmetric hydrostatic Euler equations: for 
{t,x,r) G (0,T) X T X (0,1), 

dtu^ + u^dxu^ + vfdru^ = -dxP 


(2.7) 


dxu^ +-driru^) = 0 
r 

vT\r=\ = 0 
u^\t=o = K, 


where the axisymmetric (without swirl) velocity field u := vT{t, x, r)ex+vT{t^ x, r)er 
and the scalar pressure p := p(t, x) are unknowns, Uq := Uq{x, r) is the given initial 
horizontal velocity, Cx and Cr are the unit vectors in the horizontal (i.e., x) and 
radial (i.e., r) directions respectively. 


Regarding the axisymmetric hydrostatic Euler equations (2.7) there are at least 
three fundamental problems: 

(i) the local-in-time well-posedness; 

(ii) the mathematical justification of the formal derivation; and 

(iii) the formation of singularity. 

In this paper we address all these problems. More specifically, we will first introduce 
a new sign condition, under which we will prove the local-in-time well-posedness of 


solutions to (2.7) and justify the formal derivation in the L‘‘ sense. Furthermore, 


we will also show that for a certain class of initial data, smooth solutions to (2.7) 
blow up in finite time. 

When studying the well-posedness and the mathematical justification of the 
formal derivation (i.e., problems (i) and (ii) above), one encounter the following 

Structural Difficulty: 

The radial velocity component u’" = —-dfl^{rdxU^) creates a 















6 


ROBERT M. STRAIN AND TAK KWONG WONG 


loss of one ^-derivative, so the standard energy methods typically 
fail. 

This structural difficulty can be overcome if we consider the problems under the 
following sign condition: 


( 2 . 8 ) 


1 


-dr- -d, 


1 


> (T > 0 


for some constant cr. To the best of our knowledge, the sign condition (2.8) 


IS 


new, and analogous to the local Rayleigh condition for the two dimensional flow. 
Under the sign condition (2.8)[ we can avoid the structural difficulty by eliminating 
the problematic term during the estimation. This elimination is based on the 
following nonlinear cancelation for the axisymmetric hydrostatic Euler equations 


(2.7) defining the vorticity w := drU^, for any A: = 0, 1, 2, • • • , 

f f d^u^d^uj rdrdx = ^ f ( dx{\d^u^\'^) rdrdx = 0. 

Jt Jo 2 Jq 


(2.9) 


This nonlinear cancelation is analogous to the nonlinear cancelation (2.3) stated in 
[T2] for the two-dimensional flow as well. 

Regarding the well-posedness (i.e., problem (i) above), we will solve the 
axisymmetric hydrostatic Euler equations (2.7) as long as the quantity LrU^ '■= 

-drU^ belongs to the function space 
r 

(2.10) HI ^ := -|w : T X (0,1) —>■ R; ||w||/f£ < -|-oo and tr < LrW — 

for some integer s > 4 and constant tr S (0,1), where the operator Lr := -dr, and 

r 

the norm || • ||ij£ is defined by 

(2.11) ^ \\dTLrM\L^{rdrdx)'-= [ [ 15“^L“=w|2 

11 ^ 11 ^ ^tJo 

|a|<s 

Roughly speaking, we require that the horizontal velocity component satisfies 

the sign condition 1(2.8) and the quantity LrU^ := -drU^ belongs to the function 
space 

(2.12) := {w : T X (0,1) ^ R; ||w||^/£ < +oo}. 

More precisely, we will prove the following main theorem: 

Theorem 2.2 (Well-posedness for the Axisymmetric Hydrostatic Euler Equations 


(2.7)). For any integer s > 4 and eonstant a > 0, if the given initial data Uq 


satisfies the following two hypotheses: 

(i) (regularity and sign eonditions) wq := LrU^ := -drU^ € i7£ . 


t-\ _L 

^It is worth noting that if w G then the upper bound < — can always be guaranteed 
by thb Sbbolev ,„e,b.ll.y|iBl]p,..lb=d th.l b > 0 I. ch»b„ be .biently 

there is no mathematical or physical meaning for the upper bound —. In other words, one can 

cr 

choose other upper bounds instead. We just require the boundedness for LrW so that the //£ 
norm [(2.11)| is equivalent to the weighted J/£ energy f(3. 13) | 
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(ii) (compatibility condition) / Uq rdr = A = constant, 

Jo 


then there exist a time T := T{s,a, ||wo||h£) > 0, a unique axisymmetric velocity 
field u := Cx + u^Cr, and a unique (up to an additional function oft) scalar 
pressure p such that 


(i) {u,p) solves the axisymmetric hydrostatic Euler equations (2.7) classically, 
1 


(li) w := LrU^ := J-drU^ e C{[0,T]-Hl ,,)nC\[0,T]-Hf-^), and 
(Hi) 


nl pi 

rdr = / 


Mq rdr = A. 


Furthermore, the system (2.7) is stable in the following sense: for any integer 


s' S [0,s), and any two solutions {ui,pi) := (uiCx + u[er,pi) and ( 112 ,^ 2 ) := 

1 


{u 2 ex + U 2 er,P 2 ) of system K2.7)[ if Wi := Lruf := -druf € C{[0,T]; D 


C {[0,T]; Hf ) and / uf rdr = A = constant for i = 1, 2, then there exists 

Jo 

constant Cs,s',cr,T,M > 0 such that 

(2.13) Ikl -W2||c([0,T];Hi') ^ C's.s'.^.T.mIKwi -W2)|t=ollL2(;d^d^), 


where M := max < sup ||wi||^f£, sup ||w 2 ||h£ >, and the Hf norm is defined by 

1 [O.T] [O.T] I 


( 2 . 11 ) 


Moreover, for any t € [0,T], the unique axisymmetric velocity field u{f) is also 
a vector field in three spatial dimensions. 

Remark 2.3. 

(i) (Intuition behind Norm (2.11)) Indeed, the norm (2.II) is natural for 
axisymmetric functions. First of all, viewing all axisymmetric functions as 
a special class of functions defined in a three-dimensional domain, one can see 
that the volume element is just a 2ti multiple of rdrdx because axisymmetric 
functions are independent of the azimuthal coordinate 9. Therefore, rdrdx is 
an appropriate and natural measure for the axisymmetric physical quantities. 
Regarding the measure rdr, the differential operator that can guarantee the 

fundamental theorem of calculus (B.l) is Lr '■= —dr, so it is also natural 


(ii) 


(iii) 


to measure the radial regularity for axisymmetric functions by using the 
differentiation with respect to Lr instead of dr- 

(Structure of the Proof) The proof of Theorem 12.21 is long, and will be 
separated into different parts: in Section[3]we will derive a priori estimates and 
prove the uniqueness and stability; in Section 0] we will show the existence by 
using an approximate scheme (see Subsection SH) or a reduction argument 
(see Subsection 021); finally, we will verify the C"* ^ regularity of u{t) in 
Appendix 0 

(Solvability in Non-periodic Domains) In Theorem 12.21 we solve the axisym¬ 


metric hydrostatic Euler equations (2.7) in the cc-periodic domain T x (0,1) := 
{(x,r); cc G T, r G (0,1)}. However, it is just for presentation convenience. 
For example, one may apply the result in Theorem 12.21 the finite speed of 
propagation in the cc-direction, and an argument similar in [7] to solve the 
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axisymmetric hydrostatic Euler equations (2.7) in other non-periodic domain 
like R X (0,1). For this purpose, one should also consider the uniformly local 
HI spaces: 


H 


L^uloc 


X (0,1) 


sup < OO 

lez 


where ||w||^. := ^ 

|q!|<S 

details to the interested reader. 




n-i 


rdrdx. We leave the 


Regarding the mathematical justification of the formal derivation (i.e., problem 
(ii) above), we will show the following 

Theorem 2.4 (Mathematical Justification of the Formal Derivation of the 
Axisymmetric Hydrostatic Euler Equations (2.7) I. For any e > 0, let 
and (u^jU^tP) be smooth solutions to the axisymmetric rescaled Euler equations 
(2.5) and the axisymmetric hydrostatic Euler equations (2.7) respectively. Assume 
that satisfies the sign condition (2.8) for some constant a > 0. If there exist 
constants Cq > 0 and fi € (0,4] such that 


JT JO 


+ -f |we — w|^} rdrdx 


<Coe^ 


t=o 


where w. 


:= -druf - 
r r 


and w := —drU^, then for all t € [0, T\, 




..x\2 


21 r r 12 

■ e It, — u 


-I- |we — wp} rdrdx < Ce^, 


where the constant C depends only on a, , u’’, w, Cq and T, but not on e nor 

«,<,p,). 

The proof of Theorem 12.41 will be provided in Section [SI 
Remark 2.5. 


(i) (Existence of the Axisymmetric Rescaled Euler Equations (2.5) I Since the 


axisymmetric rescaled Euler equations (2.5) is obtained directly from the 


axisymmetric Euler equations (2.3) by applying the rescaling (2.4) Therefore 


the global-in-time existence to the axisymmetric rescaled Euler equations (2.5) 


follows directly from that to the axisymmetric Euler equations (2.3) Indeed, 
the global-in-time existence of regular solutions (e.g., He := ex + U^eu € 

5 

(^([O, oo); iJ'*) n (^^([O, oo); 77®“^) for s > -) to the axisymmetric Euler 
equations (2.3) is classical; see [TB] and m for the case of as well as 


m for the case of bounded domains in for instance. 

(L^ Mathematical Justification) Since the local-in-time solution for the 


axisymmetric hydrostatic Euler equations (2.7) obtained in Theorem 12.21 
fulfills the hypotheses stated in Theorem 12.41 we justify mathematically 
rigorously in the sense that the solution obtained in Theorem 12.21 provides 
a reasonable approximation of the leading order behavior of axisymmetric and 
incompressible fluids moving in a very thin tube. 
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Regarding the formation of singularity (i.e., problem (iii) above), we will prove 
the following 

Theorem 2.6 (Finite Time Blowup for Smooth Solutions to the Axisymmetric 
Hydrostatic Euler Equations (2.7)). Consider the axisymmetric hydrostatic Euler 
equations (2.7) in the physical domain R x (0,1) := {(ai,r); a; G R and 0 < r < 1} 
instead of Tx (0,1). Let the initial data Uq :=Mg(a;,r) satisfy the following property: 
there exist a fixed position x G R and a constant horizontal velocity u® £ R such 
that 


(2.14) 


Wq(x, r) = for all r £ [0,1] 
either 92,L^Mq(x, 0) = 0 or dxLrUQ{x,l) = 0 
dxLluf^{x,r) < 0 for all r £ (0,1), 


where '■= —dr- Then there exist a finite time T > 0 such that if a solution 


(u,p) := (u^Cx + u''er,p) to the axisymmetric hydrostatic Euler equations (2.7) 
remains smooth in the time interval [0,T), then 


(2.15) lim max{||u"'(t)||Loo, ||a,^'u^(t)||i=o, ||a,^p(t)||ioo} = 00 . 

The proof of Theorem 12.61 will be provided in Section [01 
Remark 2.7. 

(i) (Periodic Domain vs. Infinite Tube) In Theorems 12.21 and 12.41 we consider 


the axisymmetric hydrostatic Euler equations (2.7) in the periodic physical 


domain T x (0,1) for simplicity. However, in Theorem 12.61 we consider the 
fluids in the infinite tube R x (0,1) instead. It is worth noting that the result 
stated in Theorem l2.6l also applies to the fluids in the periodic domain Tx (0,1) 
because one may consider the solutions in the periodic domain T x (0,1) as 
periodic flows in the infinite tube R x (0,1). 

D WdxU^WL-^.OT: WOxpWl’^. 


(ii) (Blowup) The finite time blowup (2.15) 


IS m \\u 


The hrst and last cases correspond to the infinite horizontal velocity and 
pressure gradient respectively; they are non-physical in a certain sense. The 
second case corresponds to the formation of singularity. 


3. A Priori Estimates 

The aim of this section is to derive the a priori energy estimates for proving 


the well-posedness of the axisymmetric hydrostatic Euler equations (2.7) We will 
derive these estimates in following three steps: 

(i) in order to simplify our computations, the axisymmetric hydrostatic Euler 
equations (2.7) will be rewritten as the vorticity formulation in terms of the 
new differential operator and dependent variables (3.5) in Subsection 13.11 
the weighted Hf energy estimates will be derived in Subsection 13.21 and 
an comparison principle as well as its applications to uniqueness and 
stability will be provided in Subsection 13.31 


(ii) 

(iii) 


3.1. Vorticity Formulation in New Variables. In this subsection we will 


first rewrite the axisymmetric hydrostatic Euler equations (2.7) in the vorticity 
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formulation so that we can avoid handling the scalar pressure p, which is a non¬ 
local quantity. After that, applying new dependent variables u, v, w and differential 
operator L^, we will further rewrite the vorticity formulation into a better form. 
Let us begin by introducing the vorticity formulation for the axisymmetric 


hydrostatic Euler equations (2.7) as follows. 

In order to derive estimates for system (2.7)[ one may encounter the technical 
difficulty that the scalar pressure p := p(t, x) is a non-local quantity, which and the 
horizontal velocity component are related by an integral expression: 

/•i 


dxP = -2 dx 


rdr. 


Inspired by the usual Euler equations, we can avoid handling the scalar pressure 
p by considering the vorticity formulation. For the axisymmetric hydrostatic 
Euler equations (2.7)[ the vorticitjH w is defined by w := drU^. Differentiating 
the momentum equation |(2.7)[ with respect to r and using the incompressibility 
condition (2.7)^, we obtain the vorticity equation 


Or equivalently. 


dtOJ + U^dxUJ -I- U^drOJ = - LO. 


{dt+U^dx+U^dr) (^) = 0 . 


Therefore, instead of studying equations (2.7) directly, we consider the following 


vorticity system: for (t,x,r) £ (0,T) x T x (0,1), 


(3.1) 


'll 

dtUJ + U^dxUJ + vJ'drUl = - UJ 

r 

= --drA(-) 

r \r / 

=-dxA (-) 
r \r / 

uj\t=o = Wo := drUQ, 
where the Dirichlet solver A is defined by 


(3.2) 


A{f ){t,x,r) := {\r^ - P^\-r^ - P^}f(t,x, p) pdp, 


which is the unique solution to 

'1 


--dr 

r 


drA{f)j = f for {x, r) G T X (0,1) 
A{f)\r=0,l = 0. 


Remark 3.1 (Biot-Savart Formulae). The equations (3.1)(3.1) j can be considered 


as the Biot-Savart formulae for the axisymmetric hydrostatic Euler equations (2.7) 


equations (2.7) and the vorticity system (3.1) are equivalent provided that the pole 
condition 


According to the equivalence Lemma IA.2[ the axisymmetric hydrostatic Euler 
[uat 
indi 

(3.3) 


lim ru^ {t, X, r) = 0, 

r—)-0+ 


^See Remark l2.1l for the explanation. 
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and the compatibility condition 
(3.4) 


rdr = 0 


hold. Indeed, the pole condition (3.3) is just a minor and technical assumption 
because we are looking for regular solutions. Furthermore, without loss of 


generality, we can always assume the compatibility condition (3.4) because of the 
following 


Remark 3.2 (Compatibility Condition (3.4)). Using the incompressibility condition 
(2.7) j, the boundary condition (2.7) j and the pole condition (3.3)[ we know that at 

the initial time t = 0, the integral / Uq rdr is independent of x, and hence, must 

do _ 

be a constant. Indeed, this constant is the A stated in Theorem 12.21 Therefore, 


we can always fulfill the compatibility condition (3.4) at the initial time t = 0 by 


applying a suitable Galilean transformation. Due to Lemma lA.il the compatibility 


condition (3.4) also holds for all later times t > 0 since it does initially. 


From now on until the end of Section [H we will always assume (3.3) and (3.4) 


and consider the vorticity system (3.1) instead of the axisymmetric hydrostatic 
Euler equations (2.7) Now, in order to simplify our computations, we are going to 
further rewrite the vorticity system (3.1) in terms of new differential operator and 
dependent variables as follows. 

Using 


(3.5) 


Lr ■= —d. 


and w := — = L^u, 
r 


we can rewrite the vorticity system [(3.1)| as follows: for (t, x, r) G (0, T) x T x (0,1), 

(3.6) 


' 9tW + u9a,W + vLrW = 0 

U = —LrA{w) 
V = d^A{w) 


I Wo 

w|t=o = Wo := — 
r 


where the Dirichlet solver A is defined by (3.2) In other words, Al(w) is the unique 
solution to 

J —L^Al(w) = w for {x,r) G T x (0,1) 

\ Al(w)|r= 0 T = 0. 

Indeed, there are at least three advantages for considering the new vorticity system 


(3.6) instead of the original vorticity system (3.1) 


(i) in terms of Lr, u and v, the incompressibility condition (2.7), becomes a 
simpler form: 

(3.7) dxU + LrV = 0', 


(ii) the system (3.6) is “almost” the same as the vorticity system for the two¬ 


dimensional hydrostatic Euler equations if we replace Lr by da , see Subsection 
14.21 for the details; 

(iii) the unknowns u, v and w are better quantities in the following sense: if 
u := u^ex +vLer is a smooth vector field in three spatial dimensions, then the 
horizontal velocity component is a smooth function, but the quantities vL 
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and w := drU^ may have singularities at the symmetric axis r = 0. However, 
one may check that all three quantities u, v and w are smooth even at the 
symmetric axis r = 0. As a result, we can estimate these new quantities 
relatively easier. 

Due to the above three advantages, we will study the vorticity system (3.6) instead 
of the original vorticity system (3.1) in Subsections 13.21 and 13.31 Let us end this 
subsection by providing the following three remarks: 

Remark 3.3 (Relation between u and w). It is also worth to mention that 

W = Lr-Vi 


in terms of the new differential operator and dependent variables (3.5) 


Remark 3.4 (Boundary Condition for v). Using the Biot-Savart formula (3.6)|j and 
the definition (3.2) of A, one may check that 

v|r=o,i = 0. 

Remark 3.5 (Sign Condition and Nonlinear Cancelation). Using the new differential 
operator and dependent variables (3.5) we can express the sign condition (2.8) and 
the nonlinear cancelation (2.9) as (3.11) and (3.10) in Subsection 13.21 respectively. 


3.2. Weighted Energy Estimates. The aim of this subsection is to derive the 

U) 

a priori Hf energy estimates for the quantity w := —. As a direct consequence, 

r 

we will also obtain an L°° control on as well. More precisely, we will prove 
Proposition 13.71 below. 

The main idea of the energy estimates is to estimate the weighted energy 
below instead of the original norm (2.11) Indeed, the weighted energy 


(3.13) 


(3.13) 


which is well-defined under the sign condition (2.8) (or equivalently, (3.11) 


below), is well-chosen so that it can avoid the a;-derivative loss, which is the 
structural difficulty that we mentioned in Subsection 12.21 by using the nonlinear 
cancelation (2.9) (or equivalently, (3.10) below), so the standard energy method 
works. 

These a priori estimates only reply on the fact that if (u, v, w) is a smooth 
solution to the vorticity system (3.6)[ then the quantities u, v and w will satisfy 
the following three mathematical structures: 

(i) the first order system: for {t,x,r) S (0,T) x T x (0,1), 

9(W -I- ndx'Vf + vL^w = 0 
cla-u -I- LrV = 0 


(3.8) 


v|r= 0 ,l = 0 
w|t=o = Wo := LrUo; 


(ii) the energy estimates for u and v: for any integer s > 0, there exists a constant 
Cs > 0 such that 


(3.9) 


\Dlv\\L2^rdrdx) < Cs||w||ff£ 


where := Lr := -dr and 


for any |q;| < s, 

for any |a| < s with D2 ^ 9®, 

\\hi is defined in (2.11)| and 



































AXISYMMETRIC HYDROSTATIC EULER EQUATIONS 


13 


(iii) the nonlinear cancelation: for any fc = 0, 1, 2, 


(3.10) 


dtvdtw rdrdx = 0. 


T Jo 


Remark 3.6 (Mathematical Structures). The above three mathematical structures 


can be easily verified as follows: let (u, v,w) satisfy (3.6) 


(i) Indeed, (3.6) ^ and (3.6) j are exactly the same as (3.8) ^ and (3.8) so it suffices 
to verify the incompressibility condition (3.8) j and the boundary condition 
,. However, they are just direct consequences of the Biot-Savart formulae 


(3.8) 


(3.6) 


,-(3.6)j and the definition (3.2) of A] 


(ii) Since w = L^u and dxXi + L^'v = 0, the energy estimates (3.9) just follow from 
the definition (2.11) of norm and the Poincare inequality (B.4)[_ 


(iii) Using w = Lr-u, 9a:U + Lr'v = 0, the integration by parts formula |(B. 2) and 


the boundary condition (3.8) j, we have 


rdrdx = 


T Jo 




dZoLl 1 rdrdx = 0, 


LrW := L£u > cr > 0 


since u is periodic in x. 

Under the sign condition 

(3.11) 

for some constant cr, the mathematical structures (i)-(iii) are sufficient to derive 
our a priori estimates for solutions, where 

C£“([0,r] X T X [0,1]) := {w: [0,r] x T x [0,1] ^ M; 

dldiL^r w is continuous for any non-negative integers i, j and k}. 

More precisely, we will prove the following 

Proposition 3.7 (A Priori Estimates). For any integer s > 4 and constant cr > 0, 
let w S C'([0, T]; ^), and u, v, w G ^£“([0, T] x T x [0,1]). Assume that u, v and 

w satisfy system (3.8), energy estimates (3.9)[ and nonlinear cancelation (3.10) 
Then we have the following estimates: 

(i) (Weighted Energy Estimate) there exists a constant > 0 such that 

(3.12) l|w(t)||jy, < Ijwoll^. {l + I|w(t)||^. I |lw(r)|||, dr, 

where the weighted energy || • || is defined hy 


(3.13) 




dZ 


VTyw 


L^irdrdx) 

ai^s 


and Dl := 


(ii) (L^ Estimate) there exist a time T := T{a, ||wo||//4) G (OiTj and a constant 
Ca > 0 such that for any t G [0, T], we have 
(3.14) min LrWo — CaW^oWlcit If: Lrw(T,x,r) < max -b CCT||wo||?r 4 1, 

Tx[0,l] "i- Tx[0,l] "i- 

for all (r, x,r) G [0, t] x T x [0,1]. 












































14 


ROBERT M. STRAIN AND TAK KWONG WONG 


Remark 3.8 (Equivalent Energies). As long as w G the norm (2.11) 


IS 


equivalent to the weighted energy (3.13) | becau se of the L°° boundedness of L^w. 

Therefore, the weighted energy estimate (3.12) provides a control on ||w||//£ as well. 

Proof of Proposition \S.T\ 

(i) Inequality |(3.12)| follows directly from the standard energy method because 
one can apply the nonlinear cancelation (3.10) to avoid the a;-derivative loss, 
which is caused by v = —dxLf^n. The details are as follows. 

Differentiating equation (3.8)^ with respect to Df, we have 

dtDfw + udxDfw + vLrIl£w 
=- (Z) + D^^YLrD^-'^wY 

0</3<a 

In the case that |a| < s and ai ^ s, since all D^v = —on the 
right hand side of (3.15) have no more than s x-derivatives, we can apply the 
standard energy estimates, as well as using the Sobolev inequality (B.5) and 
energy estimates (3.9)[ to obtain 

(3-16) ’^W^L^Wl^irdrilx) ^ C's||w||f^|. 


In the case that a = (s, 0), equation (3.15) becomes 
dtd^w + + vi9®LrW + d^vLrW 


(3.17) 


= -E 




dSd: 


s—fc+l_ 


-E 

k=l 


dfvdt-'^LrW. 


Now, we are going to make use of the nonlinear cancelation (3.10) 


to 


eliminate the problematic term 9®v = —dx^^Lj, ^u, which has more than s 

d: 


x-derivatives. Multiplying (3.17) by ——, and then integrating over T x (0,1) 

- LrW 

with respect to rdrdx, we have, via using the evolution equation for L^w, 


(3.18) 


1 d 

2 dt 

= -E 

k^l 

s-l 

-E 

k^l 

1 

+ 2 




LrW 


rdrdx 


dUdtw rdrdx 


d^ndtr^+^wdlw 

LrW 

Lrwdiw 


rdrdx 


LrW 
wdxW — dxULrW 


rdrdx 


iLr-wl 


|i9^wf rdrdx. 


It follows from the nonlinear cancelation (3.10) that the second integral, which 
contains the problematic term 9®v, on the left hand side of (3.18) is equal to 
0. Furthermore, using the Sobolev inequality (B.5)[ energy estimates (3.9) 


and the fact that 0 < cr < L^w < —, we know that the first two terms and 

a 

(j 

the last term on the right hand side of |(3.18)|are controlled by —-||w||?ra and 

- a ^ 
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c 

-^\\M\]iiWx^\\h{rdrdx) respectively. Therefore, we have 


(3.19) 


d 

dt 


d^w 




(rdrdx) 


< Cs,a (l + ||w||/i-*) ||w|||^3 . 


(3.20) 


Using (3.16) and (3.19) and summing up over a, we finally obtain 

^ Cs,a (l + \Mhi) ||w|||, 
since \\w\\hi < ^||w||^. 


The weighted energy estimate (3.12) follows 


directly from (3.20) 


(ii) The L°“ estimate (3.14) is a direct consequence of the weighted energy 
estimate (3.12) Indeed, inequality (3.12) implies that there exist a time 
T := T{a, ||wo||//4) £ (0,T] and a constant > 0 such that 


(3.21) 


sup ||w||^4 < C^WwqWhi- 
[O.T] 


Hence, using the evolution equation for LrW, the Sobolev inequality (B.5)[ the 
energy estimates (3.9) and the uniform bound (3.21) we have the following 
pointwise control: for any 0 < t < T, 

\dtLrw\ < \udxLrw\ + |vL^w| + Iwclj-wj + |9a;uLj.w| 


< C sup ||w||^4 

\[o.T] 


< C,, II Wo II^4, 


which implies the L°° estimate (3.14) by a direct integration. 


□ 


3.3. Comparison Principle and Its Applications to Uniqueness and 
Stability. In this subsection we will first show a comparison principle for two 
solutions. Using this comparison principle, we will also obtain the uniqueness and 


stability to the axisymmetric hydrostatic Euler equations (2.7) 

Let us begin with the comparison principle. The comparison principle is 
based on the following two properties for the difference of two solutions: 

(i) Energy Estimate: 

(3.22) ||ui ^2\\L'^[rdrdx) — f^ll^l '^2\\L^[rdrdx): 

(ii) Nonlinear Cancelation: 

nl 


(3.23) 


(vi — V 2 ) (wi — W 2 ) rdrdx = 0. 


T Jo 


Using arguments similar to part (ii) and (iii) of Remark 13.61 one can verify that 
both (3.22) and (3.23) hold for any two solutions (ui,vi,wi) and (u 2 ,V 2 ,W 2 ) of 
the vorticity system (3.6) Using energy estimate (3.22) and nonlinear cancelation 


(3.23) we will prove the following 
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Proposition 3.9 Comparison Principle). Fori = 1, 2, letwi S C([0,T];i?|)n 
C^([0, T]; i7£) such that ui, Vi and wi satisfy system (3.8), energy estimates (3.9)1 

- - - ^ - 

(3.22), and nonlinear cancelation (3.23) //w 2 satisfies 0 < a < LrW 2 < —, then 

(7 

we have the following if comparison principle: there exists a constant Ca,M,T > 0 
such that 

(3.24) ||wi — W2||c([0,T];L2(rdrd£c)) < C'ct.M.tH (wi — W 2 ) \t=o\\{rdrdx)i 


where M := max < sup ||wi||^4 , sup ||w2||//4 > . 

[lO.T] ^ [0,T] 

Proof. Let u := Ui — U 2 , v := vi — V 2 and w := wi — W 2 . Then w satisfies 
(3.25) dtW + Ull9a;W + ViLrW + Ul9xW2 + vLrW2 = 0. 

Now, we are going to eliminate the problematic term v = —dxLf^u, which causes 


a loss of cc-derivative, by using the nonlinear cancelation (3.23) Multiplying (3.25) 
w 

by - -, and then integrating over T x (0,1) with respect to rdrdx, we have 


LrW2 ' 


(3.26) 


l_d 
2 dt 

1 
2 


y/Lr^2 L'^(rdrdx) 

{dt + Midx + viLr) 
dxW2 


1 




Lr'W2 

W 


uw rdrdx — 
2 


LrW2^ 
vw rdrdx 


Iwr rdrdx 


y/LrW2 


+ M 


(rdrdx) 


L^(rdrdx)i 


where the last inequality is obtained by using the evolution equation for Lr'W 2 , 
the Sobolev inequality (B.5) energy estimates (3.9) and the nonlinear cancelation 


(3.23) 


1 


Using energy estimate (3.22) and 0 < cr < LrW 2 < —, we have 


(3.27) 


\\^\\\'^{rdrdx) ^ (rdrdx) ^ ^'<7 


cr 

W 


\/Lr-W2 


(rdrdx) 


Combining (3.26) and (3.27) we obtain 


(3.28) 


dt 


y/Lr 


(rdrdx) 


C(j,M 


y/LrW2 


Applying the Gronwall’s inequality to (3.28)[ we have 

< Ca.M.T 


\/LrW2 


(rdrdx) 


w 


y/LrW2 



(rdrdx) 


L^(rdrdx) 


□ 


which implies (3.24) because 0 < tr < LrW 2 < —■ 

- cr 

Using the comparison principle (i.e., Proposition 13.9L we can immediately 
obtain the following corollary regarding the uniqueness and stability of the 


axisymmetric hydrostatic Euler equations (2.7) 
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Corollary 3.10 (Uniqueness and Stability of the Axisymmetric Hydrostatic Euler 
Equations (2.7) I. Fori = 1, 2, let (ui,pi) := (u^ex + uler,Pi) be classical solutions 


to the axisymmetric hydrostatic Euler equations (2.7) such that / rdr = A = 

Jo ^ 

constant. For any integer s > 4 and constant a > 0, if := '■= —d^uf G 

r 

C{[0,T]; n {[0,T]; then we have the following: 

(i) (Stability) for any integer s' € [0,s), there exists a constant Cs,s',a,T,M > 0 
such that the stability estimate (2.13) holds; 


(ii) (Uniqueness) in addition, if both (iti,pi) and ( 112 , 7 * 2 ) satisfy the same initial 
data Uq, then u) = U 2 and pi = p 2 + p for some function p := p{t). 


Proof. Without loss of generality, we can assume that A = 0, otherwise, we can 
always fulfill this by applying a suitable Galilean transformation. See Remark |3.21 
for the details. 

(i) Since both (ui,pi) and {u 2 ,P 2 ) are regular enough, one may apply boundary 


condition (2.7) j, compatibility condition (3.4) and elementary argument to 


show that both solutions satisfy the pole condition (3.3) By the equivalence 


Lemma rA.2l both and {u 2 ,U 2 ,aj 2 ) satisfy the vorticity system 

where uii := druf for z = 1, 2. Equivalently, using the change of 

OJi, 


(3.1) 


variables (3.5)[ we know that (ui, v^, w^) := {uf, ru(, —) satisfies the vorticity 


system (3.6) for z = 1, 2. Therefore, we can apply the comparison principle 


(i.e.. Proposition [SHI) to obtain estimate (3.24) The stability estimate (2.13) 


follows directly from estimate (3.24) and the interpolation inequality (B.7) 


(ii) Since both (zzi,pi) and {u 2 ,P 2 ) satisfy the same initial data Uq, one can check 
that wilt^o = W 2 |t= 0 ) where w^ := -druf for z = 1, 2. Thus, applying the 
stability estimate (2.13)[ we obtain wi = W 2 . Using the Biot-Savart formulae 


we also have ui = U 2 and vi = V 2 . In terms of the original 


variables, we obtain {ui,u[,u!i) = ( 712 ,^ 21 ^ 2)1 and hence, by the equivalence 
Lemma [A.21 pi = p 2 + p for some function p = p{t). 


□ 


4. Existence 


The aim of this section is to provide two independent constructions of the 

In Subsection 


solutions to the axisymmetric hydrostatic Euler equations (2.7) 


10 we will first introduce an approximate scheme that keeps all a priori estimates 
derived in Section [31 and then construct the solution as the limit of approximate 
systems. In Subsection 14.21 we will construct the solution by a reduction argument. 


4.1. Existence via Approximate Scheme. In this subsection we will first 
introduce an approximate scheme. Using this approximate scheme, we will outline 


the proof of existence to the axisymmetric hydrostatic Euler equations (2.7) 
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Let US begin by introducing the approximate systems as follows: for any positive 
integer TV, we consider, for (t,x,r) G (0,T) x T x (0,1), 


(4.1) 


' dtWjsf + = 0 

ujv = —PNLrAiyff^) 
VAT = PNdxA{wN) 


I Wo 

WJV t=o = Wo := —. 

r 


Here, the Dirichlet solver A is defined by (3.2) and the projection operator Pjv is 
defined by 


PNf{t,x,r) 

(4.2) ^ 

:= aoH,rj + 

k^l 

where the coefficients are given by 


N N 

:= ao(t,r) + E ak(t,r) ^•\/2cos2fc7rx^ +y^6fc(t,r) ^•\/2sin2fc7ra;^ , 


fc=i 


ao{t,r) := / f{t,x,r) dx 

Jt 

afc(t, r):= J f {t,x,r) cos 2kTTx^ dx for all fc = 1, 2, • • • 
bk[t,r) := J f{t,x,r)(^-\/2sm2kTTx^ dx for all fc = 1, 2, • • •. 


The main advantage of the approximate system (4.1) is that for any fixed TV, 


the system (4.1) does not have the loss of x-derivative (i.e., the structural difficulty 
mentioned in Subsection l2.2l) because the projection operator Pjv regularizes v^r in 
the x-direction. As a result, for any fixed TV, one may construct the unique local-in¬ 


time solution to the approximate system (4.1) by classical methods. For example, 
one may first construct an iterative sequence of linearized solutions to (4.1) and 


then prove the convergence of these linearized solutions by using the energy method. 
Since the details are standard, we leave this to the interested reader. 

It is worth noting that if we only apply the standard energy method to solve the 


approximate system (4.1) then the life-spans of the solutions may depend on TV. 


This is due to the fact that the regularization effect of the projection Pjv becomes 
weaker and weaker for a larger and larger TV. Therefore, in order to solve the 


approximate system (4.1) in a uniform (in TV) life-span, one must derive a priori 
energy estimates without using the regularization effect of Pn- This can be done 
under the sign condition (3.11) More precisely, one may check that the solution 
(uat. Vat, Wat) of (4.1) indeed satisfies the mathematical structures (3.8)|(3.f0)| and 
hence, we can also apply the Proposition 13.71 to the solution (uat, vat,wat). As 


a result, weighted energy estimates (3.12) and L°° estimate (3.14) also hold for 


Wat as well. Since estimates (3.12) and (3.14) are uniform in TV, one may apply 
the standard continuous induction argument to show that for any integer s > 4, 
TV = 1, 2, • • • and ct > 0, if wq G 2 o-j there exist a uniform (in TV) life-span 

T > 0 and a solution wjy G (^([O, T]; „,) fl (^^([O, T]; P£“^) to the approximate 
system (4.1) Furthermore, the solutions wjy are uniformly bounded: there exists 
a constant M > 0, depending on s, a, T and ||wo||Ar£ but not on TV nor wat, such 
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that 


(4.3) 


l|wAr||c([0,T];H|) < M. 


Finally, in order to solve the vorticity system (3.6) we have to prove the 


convergence of w^r and the consistency of the limit. Regarding the convergence, one 
may follow a similar argument as in Subsection 13.31 to verify the convergence. 
However, for any two different approximate solutions, the nonlinear cancelation 
(3.23) does not hold in general. Indeed, we only have the almost nonlinear 
cancelation: for any integers iVi > N 2 , 


(4.4) 


If 


{vNi - vatJ {wni - watJ rdrdx 


— ^2 IKwi '''^^2llL2(rdrda:) : 


where the constant Cm depends on the constant M mentioned in inequality (4.3)[ 
but not on and N 2 . As a result, if we follow the argument in the proof of 
Proposition 13.91 and use the fact that w^vJ^^q = then we can show, for 

any integers Ni > N 2 , 


IKaTi - ^N2\\c{[0,T]-L'^{rdrdx)) < 


Ca^M,T 

Ni ■ 


This proves the C ([0, T]; convergence. By the interpolation inequality (B.7) 


we 


obtain the C ^[0,T]; ^ convergence of the approximate solutions w^r as N goes to 

00 for any integer s' G [0, s). Since we have the strong convergence, the consistency 
is obvious. Lastly, using the standard regularizing initial data argument, one may 
also construct w G C ([0, T]; iL£). Since the estimate (4.3) also holds for w, we 
can verify by using the evolution equation for that if wq G 20 -j then w G 
C (\0,T];HI for some T > 0. This and the equivalence Lemma [A. 2 1 show the 
existence stated in Theorem 12.21 

Let us end this subsection by verifying the almost nonlinear cancelation (4.4) as 
follows: 


Proof of Inequality (4.4) 


First of all, we can rewrite 


(4.5) 


where 


{vNi - vats ) (watj - WN 2 ) rdrdx 


T Jo 


h = 


h = 


/ / (vAi - ^Ai5a:A(wAfJ) (waTi 

Jf Jo 

/ / (TVi9a;A(WA2) - VA2) (WAi 

Jt Jo 


/ 1 +/ 2 , 


watj) rdrdx 
WA 2 ) rdrdx. 
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Using VNi = PNidxA{wNi) and Wjy^ — wjVa = —L^A^wni —WN 2 ), we can apply 
the integration by parts formula (B.2) to the integral Ii, and obtain 


h = - / PNidxA{wNi -wn 2 )LIA{wni - wa?J rdrdx 

Jt Jo 

= PNidxLrA{w]sr^ - wjsr^) LrA{wN^ - wjsr^) rdrdx 

(4-6) Jo 


dx \PNiLrA{wNi - watJI^ rdrdx 


T Jo 


= 0 . 


On the other hand, using iVi > N 2 and vn 2 = PN 2 dxA{wN 2 )^ we can estimate I 2 
as follows: 

^2 = / / {{Pni - PN2) dxA{wN2)} (watj - watJ rdrdx 
Jt Jo 

^ II (Pni - PN2)dxA{wN2)\\L'^(rdrdx)\\^Ni “ Wat^ || (rdrdx) 


(4.7) 


C 3 

- ]^ll - PN 2 ) i9a:-4(wAr2)||L2( rdrdx) ||L^(rtir(ia:) 

c 

— jyT N2^L'^(rdrdx) 

^ ^2 N2^L‘^(rdrdx)-> 


where the constant Cm depends on the constant M mentioned in inequality (4.3) 
Combining (4.5)|(4.7) we prove inequality |(4.4)' 


□ 


4.2. Existence via Reduction. The aim of this subsection is to introduce a 


reduction argument which reduces the vorticity system (3.6) to be the vorticity 


system for the two-dimensional hydrostatic Euler equations. As a result, we can 


obtain the solutions to the axisymmetric hydrostatic Euler equations (2.7) by using 
that to the two-dimensional hydrostatic Euler equations, whose well-posedness 
were shown in m previously. 

The main idea of the reduction argument is to use the cross-sectional area a 
instead of the radius r as an independent variable. To be more precise, let us begin 
with the following change of variable: define 


(4.8) 


1 2 
a := —r 
2 


which is equivalent (up to a factor of 27r) to the area of a two-dimensional disk of 


radius r. It follows from (4.8) that 


da = —dr ='■ Lr and da = rdr. 
r 


(4.9) 
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Thus, using the change of variable (4.8) we can rewrite the vorticity system (3.6) 
as the following system: for (t,x, a) £ (0,T) x T x (0, -), 

9tw + u9xW + vda'w = 0 

u = -i9a^(w) 

V = dxA{w) 
w|t=o = Wo, 

where the Dirichlet solver A is given by 

1 

1 


(4.10) 


(4.11) 


A{f){t,x,a) := -- 


|a — a| — 0 — 0 + 4oo}/(t, x, a) da. 


In other words, A{f) defined in (4.11) is the unique solution to 

-dlAif) = / in T X (0, i) 

Rescaling the interval (0, -) to be (0,1) and viewing the variable a as the length 

instead of area, one can recognize that the vorticity system (4.10) is exactly the 
same as the vorticity systeir0 for the two-dimensional hydrostatic Euler equations. 
Furthermore, by (4.9)[ the sign condition (3.11) can be written as 

daW = dlu > cr > 0, 

which can also be seen as the local Rayleigh conditioifl for the two-dimensional 
hydrostatic Euler equations. 

As a result, it follows from Theorem 2.5 in [12] that for any integer s > 4 
and constant a G initial data wg G ^ 

vorticity system (4.10) has a unique solution w G C{[0,T]-, H^{T x (0, i))) fl 
(^^([O, T]; X (0,-))), where the life-span T depends on ||wo||h«, s and a 

only. Here, the function space x (Oj ^)) is defined by 


1 . 


1, 


HliJ X (0, -)) := w G H^{T x (0, -)); 0 < a < ^.w < 


1 


and 77®(T x (0, -)) is the standard space defined on T x (0, -) equipped with 
the standard H‘^ norm: 


(4.12) 




-E 

lal^s 


dxda. 


Using (4.9)1 one can verify that the standard norm (4.12) with respect to x 
and a is equivalent to the Hf^ norm defined in (2.11) and hence, the function 

U. . 1 


spaces H^{T x (0, -)) and H^(T x (0, -)) are equivalent to HI defined in (2.12) 


"^See system (2.1) in m for instance. 
^See inequality (2.2) in m for instance. 
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and HI ^ defined in (2.10) respectively. Therefore, the existence stated in Theorem 
12.21 follows directly from Theorem 2.5 in m and the equivalence Lemma IA.2I 

5. Mathematical Justification of the Formal Derivation 

The aim of this section is to provide a rigorous mathematical justification of the 
formal derivation of the hydrostatic Euler equations (2.7) In other words, we will 
verify the formal derivation introduced in Subsection 12.11 by proving Theorem 12.41 
Let us begin by rewriting the axisymmetric rescaled Euler equations (2.5) in 

terms of the differential operator L,. := —dr and new dependent variables u^, v^ and 

r 


Wg, which are analogous to the variables u, v and w introduced in (3.5) previously. 
Define 

£2 

(5-1) Ue := Mg, Vg := rUg, and Wg := - ^d^Ve- 


Using (5.1) we can express the axisymmetric rescaled Euler equations (2.5) as the 


following equations: for (t,x,r) £ (0,T) x T x (0,1), 

5tUg + Ug9a;Ug + VgL^Ug = -dxPe 

+ nA + VgLr I = -drPe 
dx^e + LrVe = 0 
Ve|r=0,l = 0 
Ug|t=o = UgO := Mgo- 


(5.2) 


In terms of the change of variables (5.1) the corresponding vorticity system, called 


the axisymmetric rescaled vorticity system, becomes: for (t, x, r) £ (0, T) xTx (0,1), 

9tWg + Ug9xWg + VgLj.Wg = 0 

9a; Ug + LrVe = 0 

' v.|„,.=0 

Wg|i=o = WgO := LrUeO- 

Now, the question becomes as follows: for any solution (u,p) := (u^ex + u’"er,p) 
to the axisymmetric hydrostatic Euler equations (2.7) does there exist a sequence 
{(ug,Vg, Wg)}g>o such that (ug,Vg,Wg) converges to (u, v,w) := (u^,ru^, LrU^) as 
e goes to 0“''? Under the sign condition (3.11)[ the answer is affirmative in the 
following sense: 

Proposition 5.1 (Mathematical Justification of the Formal Derivation). Under 
the hypotheses of Theorem \2.4[ if there exist constants Cq > 0 and j3 £ (0,4] such 
that 

„ „i f 

Vg — V 



Ug - up + 


+ Wg - w ■ 


rdrdx 


< 


t=o 


/T JO 

then for all t £ [0, T], 
(5.4) 


where the constant C depends only on a, u, v, w, Co and T, but not on e nor 

(Ug, Ve,Wg). 


f 


2 'j 

|ug-up+e2 

Vg — V 

+ Wg - wp 

r 

1 

i 
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It is worth noting that Theorem 12.41 is a direct consequence of Proposition 15.11 
Following the approach in [2] , we will prove Proposition 15.11 by using the entropy 
method. Under appropriate modifications, our proof is similar to the proof in [5], 
in which the formal derivation of the two-dimensional hydrostatic Euler equations 
was rigorously justified. These modifications are necessary because the structure 
of the axisymmetric rescaled Euler equations (5.2) is differen10 than that of the 
two-dimensional rescaled Euler equation^^. 


Proof of Proposition 15.11 In this proof we will first define a convex functional that 
is equivalent to the left hand side of (5.4) Then we will derive a growth rate 
control on this convex functional. The main difficulty of this proof is that the 
terms involving v^ — v may cause a loss of one e, which is corresponding to the x- 
derivative loss caused by the hydrostatic limit. However, this e loss can be avoided 
by a well-chosen convex functional due to the nonlinear cancelation. 

First of all, we are going to define the convex functional by using the following 
convex (in w) function: 

Lemma 5.2 (Convex Function). For any given u and w satisfying the sign 
condition (3.11)[ there exist a constant k, and a smoot/0 and strongly conve^ (in 
w ) function E : [0, T] x T x R —> R such that for all (t, x, r) G [0, T] x T x [0,1], 


(5.5) 


dlF {t,x,w{t,x,r)) = 


u(t, x,r) — K 
Lrw(t, x, r) 


> 1 . 


Ffere, the norm of F depends on u, w, and the a stated in the sign condition 


(3.11), but the constant k depends on u and w only. 


It is worth noting that both k and F stated in Lemma 15.21 are independent of e 


and (ue,Ve,We). Furthermore, the constant 1 stated in (5.5) is not a crucial value: 


the proof of Proposition 15.11 will also work if it is replaced by any other positive 
constant. Assuming Lemma l5.21 which will be shown at the end of this section, for 
the moment, we can define the convex functional by using the convex function F 
as follows. 

For any given (u, v,w), we define the following convex functional 
(5.6) L,{t):=LkAt)+LcAt), 

where the kinetic energy part and the convex part Lc,e are given by 


(5.7) 


Lk,e{t) := 



Ue - Up + 


Vg — V 

n 

r 

1 


rdrdx 


Lc,c(t) '=1 [ [ {F{t,x,w^) - F{t,x,w) 

^ JT Jo 

— dyxFlt, x, w)(we — w)} rdrdx. 


®More precisely, the axisymmetric rescaled Euler equations | (5.2) | is not equivalent to the two- 
dimensional rescaled Euler equations even if we apply the reduction argument stated in Subsection 

[U 


^See equations (9.2) in m for instance. 

^Here, “smooth” means that for any fixed {t,x) £ [0, T] x T, the functions F, dtF and dxF 
are twice continuously differentiable with respect to w, and the partial derivatives dtd^F, 

dxd-wF and dxd^F are bounded in [0,T] x T x M. 

^That is, there exists a constant c > 0 such that d^F > c. 
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Using the smoothness and convexity of F, one may verify that the convex 
functional Le is equivalent to the left hand side of (5.4)[ so in order to prove 
Proposition 15.11 it suffices to show 

Leit) < 

for /3 S (0,4] provided that L^{0) does initially. 

A direct computation, which we leave to the interested reader, yields 


Lemma 5.3. 

(5.8) 

where 

h ■■= -- 
2 

F := 


dt 

dt 


Lk,e — Ii + F F Y + R 


Lc,t — F F F F F F X, 


da;U < lue - up + 


T Jo 


Ve — V 

n 

r 

J 


rdrdx. 


IL (?) <"•-">(■ 




rdrdx 

2 


rdrdx, 


F '■= — 


dxwd^F{t, X, w) (ue — u) (wg — w) rdrdx, 


T Jo 


h ■= [ [ {dtF{t, X, We) - dtF{t, X, w) - dtd^F{t, x, w) (w^ - w)} rdrdx, 

Jt Jo 

■ = ne{dxF{t,x,wF - dxF{t,x,w) - dxd,„F{t,x,w) {w^ - w)} rdrdx, 

JT JO 


F ■■= 

X ■=- 

Y := 


d^F(t, x, w)LrW (ve — v) (wj — w) rdrdx. 


T Jo 


/T Jo 
,2 


u (ve — v) (wg — w) rdrdx, 
1 


R:——e‘^J j (clt+vLr)^—^ 


rdrdx. 


Applying Lemma 15.31 we have 
(5.9) 


dt dt 

where Z := X FY, that is, 

nl 


d d d T rr r, 

= ^FFZfR, 


i=l 


Z = - 


T Jo 


{9pF(t, X, w)LrW — u} (ve — v) (wg — w) rdrdx. 


It is easy to check that for alH = 1, 2, • • •, 5, 

(5.10) \F\ < L„ 

so we only have to control Z and R. Indeed, both Z and R have the factor Ve — v, 
which is problematic because it may create a loss of one e. Therefore, in order to 
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obtain good estimates for Z and i?, one should eliminate the problematic factor 
Ve — V. This can be done by making use of the well-chosen convex function F and 
the integration by parts argument as follows. 


First of all, to deal with Z, we apply (5.5) to rewrite 


Z = K 


(ve — v) (we — w) rdrdx. 


T Jo 


Using the identity 


We - W = Lr{\Xe - u)- ndx^e 


the integration by parts formula (B.2) and the incompressibility conditions (3.7) 
and (5.3),, we have 


(5.11) 


Z = K / (ve — v) Lj. (ue — u) rdrdx 

Jt Jo 

I I ( dx ( 1 rdrdx 


T Jo 


/iM)) 


dx < lue - u| - 


rdrdx 


n Jo 
,2 


> 

1 

> 

n 

r 

i 


rdrdx 


= —Ke 


/iMC) 
/i'M)) 


r 

Vf — V 


rdrdx 

rdrdx 


since all Ug, Ve, u and v are periodic in x. It is worth noting that in (5.11) 


we 


applied the nonlinear cancelation argument to eliminate the leading order term. 
Combining (5.9)|(5.lT) we have 


(5.12) 


dt 


Le — ^ ( li ~\- Z -\- R ^ Co-, \1,V,W Le + 


i=l 


where R is defined by 


R := + Kdx + vL, 

JT JO 


■>( 7 ) 


rdrdx. 


Finally, we can estimate i? by a simple integration by parts argument. Define 


pT 

q{t,x,r):= J {dt + ndx + vLr) 


dr. 


{t,x,r) = {t,x,r) 
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Then using the integration by parts formula (B.2)[ and the incompressibility 
conditions (3.7) and (5.3),, we have 


(5.13) 


R = —e^ / LrQ (ve — v) rdrdx = / / qLr (v^ — v) rdrdx 

Jt Jo Jt Jo 

qdx (ue — u) rdrdx = / / dxQ (u^ — u) rdrdx 

Jt Jo 


T Jo 
< CgL, + e^. 


Combining (5.12) and (5.13) we obtain 

d 


dt 


Tg ^ Cgr u,v,w Tg C . 


Applying the Gronwall’s inequality, we have 

Le{t) < {Tg(0) + eH} *, 


which implies (5.4) for /3 G (0,4]. 


In order to complete the proof of Proposition EH we will show Lemma 
follows: 


□ 


as 


Proof of Lemma 1 5. ill Following a similar argument in [5] , we will construct the 

convex function F as follows. Since —drW = L^w > ct > 0, we know that for any 

r 

fixed (t,a;), the mapping r i—>■ w{t,x,r) is strictly increasing. Therefore, let 
n := {(t,a:,w) G [0,T] x T x R; w{t,x,0) < w < w{t,x, 1)}, 
and we can define a function i? : 11 —>■ R such that 

R{t,x,w):=r if w{t,x,r)=w. 

In other words, for any {t,x,r) G [0,T] x T x [0,1] and {t,x,w) G 11, we have 
r = R{t,x,w(t,x,r)) and w = w(t,x, R{t,x,w)). 

Now, for any fixed constant k G R, we can define, via using this mapping R, 

u(t, X, R{t, X, w)) — K 


G{t, X, w) := 


LrW{t, X, R{t, X, w)) 


for all {t, X, w) G H. We can further define F by integrating G twice with respect to 
w, and extending F smoothly in the variable w. The extension here is not unique 
and not important in the proof of Proposition 15.11 To obtain inequality (5.5)[ k 
must be chosen appropriately, for example, we can choose k := infu— supL^w. 
Here, the choice of k is not unique as well. □ 


6. Formation of Singularity 

In this section we will prove the blowup result (i.e., Theorem [SHI for the 


axisymmetric hydrostatic Euler equations (2.7) The proof is based on the reduction 


argument introduced in Subsection 14.21 and the recent blowup result in m- 
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Proof of Theorem \2.(A Using the new differential operator and dependent variables 


(3.5) we can rewrite the axisymmetric hydrostatic Euler equations (2.7) as follows: 
for {t, X, r) G (0, T) X K X (0,1), 

9tu + u9a;U + vLf.u = —dxP 
dxn + LrV = 0 
^ v|,=o.i=0 

u|t=o = Uo. 

Similar to Subsection 14.21 we can apply the change of variables 

a := r^ and v := 2v 


to rewrite (6.1) as the following equations: for {t,x,a) G (0,T) x R x (0,1), 

dtVL + u9a;U + vi9au = -dxP 

dxU + daV = 0 

V|a=0.1 = 0 
U|t=0 = Uo, 


( 6 . 2 ) 


which is exactly the same as the two-dimensional hydrostatic Euler equations if we 
see a as the length instead of area. 

Similarly, under the same change of variables, one may check that the hypothesis 


(2.14) becomes 


(6.3) 


uo(i, a) = u for all a G [0,1] 
either cl2;9aUo(a:, 0) = 0 or cl2;9aUo(i, 1) = 0 
dxd ? uo { x , a ) < 0 for all a G (0,1), 

where u := is the given constant horizontal velocity. Without loss of generality. 


we can further assume that (6.3), is just 

(6.4) 9a;i95Uo(x,0) = 0; 

otherwise, we can apply the change of variables a := 1 — a and v := — v as well. 

Since the initial data uq satisfies the hypotheses (6.3) p (6.3)1^ and (6.4)[ by 
Theorem 2.1 in m, there exists a finite time T > 0 such that if a solution (u, v,p) 
to the two-dimensional hydrostatic Euler equations (6.2) remains smooth in the 
time interval [0,T), then 

(6.5) lim max{||u(<)||Lcx>, |ia,,u(t)||L~, ||9a,p(t)||L=o} = oo. 

t-s-T- 


Indeed, (6.5) is equivalent to (2.15) under the above change of variables, so Theorem 
is just a direct consequence of Theorem 2.1 in [HI- □ 


Appendix A. Basic Properties for the Axisymmetric Hydrostatic 

Euler Equations 

In this appendix we will study the basic properties for the solutions of the 


axisymmetric hydrostatic Euler equations (2.7) More precisely, we will show that 


the average horizontal velocity is conserved, and the axisymmetric hydrostatic Euler 
equations (2.7) and its vorticity system (3.1) are equivalent. 
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Lemma A.l (Conservation of the Average Horizontal Velocity). Let {u,p) := 
{u^Sx + u^Brjp) be a classical solution to the axisymmetric hydrostatic Euler 
equations (2.7) in the periodic physical domain Tx (0,1) such that the pole condition 


(3.3) holds. Then 


(A.1) 


[ u^rdr= [ 
/o Jo 


Uq rdr. 


Proof. Integrating (2.7)| and (2.7) j with respect to rdr over (0,1), we have, after 
applying integration by parts, the boundary condition (2.7)1, and the pole condition 


(3.3) 


(A.2) 


% [ rdr + dx [ 

Jo Jo 


u-f rdr = -^dxP 


dx J rdr = 0. 


It follows directly from (A.2), that / rdr is independent of a;. Thus, integrating 

_ Jo 


(A. 2)1 with respect to x over T and using the x-periodicity of the functions and 
p, we obtain 


dt 


rdr = 0, 


which implies the identity (A.l) 


□ 


Lemma A.2 (Equivalence Lemma). We have the following: 

(i) Let be a classical solution to the vorticity system (3.1) Define 

(A.3) 


p-.= -2 


rdr + p, 


where p := p[f) is any arbitrary function oft. Then {u,p) := {u^Cx + u''er,p) 
solves the axisymmetric hydrostatic Euler equations (2.7)[ and satisfies the 


pole condition (3.3) and the compatibility condition (3.4) 


(ii) Let the initial horizontal velocity Uq satisfy the compatibility condition (3.4) 
and {u,p) := {u^Cx + u^er,p) be a classical solution to the axisymmetric 


hydrostatic Euler equations (2.7) If satisfies the pole condition (3.3) and 
uj := drU^ € C'^([0,T] X T X [0,1]), then (u^,u’',uj) solves the vorticity system 


(3.1)[ and satisfies the compatibility condition (3.4) for all later times t > 0. 


Proof. The proof is just a direct checking and will be outlined as follows. 


(i) It follows directly from the Biot-Savart formulae (3.1)(3.1), and the 
of A that u) = drU^, and (d, p) satisfies the incompressibility 


(3.2) 


(2.7) 


definition 
condition 

well as the compatibility condition (3.4) It remains to show that (d,p) also 


the boundary condition (2.7)|,, the pole condition (3.3) 


as 


I- 


satisfies the momentum equation (2.7), and the initial condition (2.7), 

Substituting w = into the vorticity equation (3.1)| and using the 

incompressibility condition (2.7),, we have 


dr {dtu^ + U^dxU^ + U^'drU^} = 0 , 


(A.4) 


and hence, a direct integration yields, via using the boundary condition (2.7) 
dtu^ + U^dxU^ + U^drU^ = q, 
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where q := {dtu^ + u^dxU^} |r=i is a function of t and x only. To see that q 
is actually equal to —dxP, we integrate (A.4) with respect to rdr over (0,1), 
and obtain 


(A.5) 


-q= u^dxU^ + u'"drU^ rdr = d. 


f 

Jo 


rdr = -^dxP- 


Here, we applied (3.4) in the first equality. In the second equality, we 

The third equality 


applied integration by parts, (2.7)(2.7), and 
follows directly from formula (A.3) Combining 


the momentum equation (2.7),^. 


HL 

(A.4) and (A.5)[ we verify 


Finally, using the Biot-Savart formula (3.1) j and w = drU^, we have 


it=o 


= --drA(—) = --drA 

r \ r / r 


-drU^ = Ug, 


where the last equality holds under the compatibility condition (3.4) This 


verifies the initial condition (2.7) 


(ii) First of all, a direct differentiation of initial condition (2.7) 


and completes the proof of part (i). 

with respect 


to r yields the initial condition (3.1) j. Next, differentiating the momentum 


equation (2.7)| with respect to r and using the incompressibility condition 


(2.7),, we obtain the vorticity equation (3.1) 


According to Lemma lA.ll the compatibility condition (3.4) holds for all 
later times t > 0 because it does initially. Finally, the Biot-Savart formulae 
(3.1)(3.1) j can be verified by using the standard Green’s function technique 
because and vJ' are the unique solutions to 


- {ru^) = - 



1 = 0,1 = 0 


respectively. 


□ 


Appendix B. Properties for the Operator Lr 


The aim of this appendix is to provide basic properties for the operator Lr 
and its corresponding measure rdr. In particular, we will discuss the fundamental 
theorem of calculus, the integration by parts formula, the Poincare type inequality, 
the Sobolev type inequality and the interpolation inequality. 

Let us begin by stating without proof the following calculus facts for Lr and rdr: 

Proposition B.l (Calculus Facts for Lr and rdr). Define the operator Lr := —dr- 

r 

Then we have 

(i) (Fundamental Theorem of Calculus) for any continuously differentiable func¬ 
tion f, for any a, b G [0,1], 


Lrf rdr; 


(B.l) 


f{b) = f{a) -b 
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(ii) (Integration by Parts Formula) for any continuously differentiable functions 
f and g, for any a, b € [0,1], 


(B.2) 


pb pb 

/ fLrg rdr = [fg]l^^ - / gUf rdr. 

J a J a 


The proof of Proposition IB.ll is elementary, so we leave this to the reader. As 
long as an operator L and its corresponding measure /i satisfy the fundamental 
theorem of calculus, the following Poincare type inequality follows immediately. In 
particular, inequality (B.4) below holds when L := and dg, := rdr. 

Proposition B.2 (Poincare Type Inequality). Let L be an operator and g be a 
measure on [0,1]. Assume that the operator L and its corresponding measure g 
satisfy the fundamental theorem of Calculus: for any function f in the domain of 
L, for any a, b € [0,1], 

(B.3) 


m = f{a) + f Lf dg. 
J a 


If there exists a point tq G [0,1] such that f^r^) = 0, then the following Poincare 
type inequality holds: 


(B.4) 


\f\Ug<g{[0,l]f \Lf\Ug. 


Proof. It follows from /(tq) = 0 and (B.3) that for any r G [0,1], 


fir) = f 

J rn 


Lf dg, 


and hence, by the Cauchy-Schwarz inequality, 


l/WP< 


\Lf\ dg 


<m([0, 1])/ \Lf? dg. 
Jo 


Integrating the above inequality with respect to the measure g over [0,1], we prove 
inequality (B.4)[ □ 


Furthermore, we also have the following Sobolev type and interpolation inequal¬ 
ities: 

Proposition B.3 (Sobolev Type Inequality). Define Lr := —dr. There exists a 

r 

universal constant C > 0 such that 
(B.5) ll/ll 

L°^{rdrdx) — ^ { ll/llL^(rdrrfa::) “b II (rrfrtia;) “b IIII (rdrrfa::) } : 

for any function / : T x (0,1) —>■ R. 

Outline of the Proof. Let a := —r^. Then da = —dr = Lr and dadx = rdrdx. 

_ 2 r 

Therefore, inequality (B.5)| is equivalent to 

(B-6) ||/||L“(dada:) < C {\\f \\l:‘{ dadx) + \\dxf\\L:‘(dadx) + \\dlf\\L^(^dadx)} , 

which can be shown by elementary methods. For instance, see Lemma B.2 in |13] 


for the proof of (B.6) in a similar domain. 


□ 
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Proposition B.4 (Interpolation Inequality). For any integers s > 0 and s' G [0, s), 
there exists a constant Cs s' > 0 such that 


(B.7) 


\Hf — 


ll.^llL2(rdrda:)’ 


where the norm || • Whs is defined in (2.II) 


Outline of the Proof. Similar to the proof of Proposition IB.31 we define a := 2 ^^- 

Since da = —dr = Lr and dadx = rdrdx, we have 
r 


Ih?. 


■=Y. ! f rdrdx 

I 1^ J'fJQ 


|a|<s 


= E / / \dTdTf? dadx =: \\fWl^(riad.y 

I 1^ JTJO 
|ck|<s 

where || • Wn^'^dadx) is the standard H" norm with respect to x and a. As a result, 
the interpolation inequality (B.7) is equivalent to 


\H‘'{dadx) ^ ^e,s'\\f\\H'^{dadx)\\f\\L2{dadx}’ 


which is just the the standard interpolation inequality. Thus, (B.7) holds. 


□ 


Appendix C. Properties for Three Dimensional Smooth Vector Fields 

The aim of this appendix is to provide a brief explanation why the unique solution 
obtained in Theorem 12.21 is a smooth vector field in three spatial dimensions. We 
will first quote the characterization result in |10] , and then, apply this to verify the 
regularity of our solution. 

In [lOj , Liu and Wang gave a characterization of smoothness of axisymmetric and 


divergence-free vector fields in terms of the pole conditions (see (C.l), below), which 


are compatibility conditions at the axis of symmetry in a certain sense. Using the 
notation in this paper, we can restat43 a special case of Lemma 2 of m regarding 
axisymmetric vector fields without swirl as follows: 


Proposition C.l (Special Case of Lemma 2 in [TO])- Denote r := \/ y^ F . Let Cx 
and Cr be the unit vectors in the horizontal (i.e., x) and radial (i.e., r) directions 
respectively. As a three-dimensional vector field, u := u'''{x,r)ex u"'{x,r)er is 

divergence-free and belongs to C^(]R^) if and only if there exists a scalar-valued 
function if := fi{x,r) G C^“''^(]R x K+) such that 


(C.l) 


u"' = Lr{rif), u'’ = -dxtp, 
= 0 for all 0 < 2m < k, 


where Lr '.= -dr 
r 


^*^The way that we state Proposition [CT] is slightly different than the original way in HOI . For 
the original statement, please refer to Co]. 
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Now, we are going to apply ProDOsition lC.il to provide a better understanding 


of our solutions to the axisymmetric hydrostatic Euler equations (2.7) Up to 
a Galilean transformation (c.f. Remark 13. 2 L the axisymmetric velocity held u 
obtained in Theorem 12.21 belongs to 


= 0 , 


:= := w'^(a;,r)ea; + M’'(a;,r')er; J rdr = 0 , 


dxU^ + L^{ru'') = 0, and w := G >, 


1 


where Lr := —dr, and the function space HI is dehned in (2.12) Indeed, as a 

three-dimensional vector held, our solution u is C‘‘~^ at the axis of symmetry (i.e., 
r = 0) as well. More precisely, we have the following 


Proposition C.2 (Regularity of Vector Fields). Let r := \Jand LI := 
{{x,y,z)-, cc G T, y^ -\- z^ < l). For any integer s > 2, 

(C.2) VLCC^-2(U). 

Outline of the Proof. Dehne (j){x,r) ■= — u^{x,f) rdf. Then we have 

J r 

= Lr4>, ru^ = -da:4>, and <(>1^=0,i = 0- 

In other words, — is equal to the scalar-valued function stated in Proposition 
r 


1C. II Furthermore, using the Sobolev inequality (B.5) the Poincare inequality (B.4) 
and the standard density argument, one may verify that ft is C^~ outside any 

neighborhood of r = 0, so is "0 = —. 


According to ProDOsition lC.il it remains to verify the pole condition (C.l),. It 
follows from a direct computation and L’Hopital’s rule that 


lim = lim 9^ ( ^ ) = lim V ^ 

^ 0 + r^ 0 + \r J r^ 0 + r-?+l ^ 

i —0 


ll 


= ito -tlULa,, (^ (zir^ 


1 


r^0+ (j -I- I)rl 


\i =0 


j + 1 r^0+ 


lim 


Therefore, it sufhces to show that lim = 0 for all non-negative integer 

]—> 0 + 

m < 2 ~ Since dr = rLr, a direct computation yields 


(C.3) 


2m+l 

E ' 




for some constants am,i- Using the Sobolev inequality (B.5) and the Poincare 
inequality (B.4)[ o ne may check that the absolute values of all on the right 
hand side of (C.3) are bounded by ||w||//?, and hence, passing to the limit r —>■ 0^ 


(C.3) we obtain lim = 0 for all non-negative integer m < -1. This 

0 + _ 2 

□ 


completes the proof and verifies (C.2) 
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